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We observed atmospheric gamma-rays around 10 GeV at balloon altitudes (15~25 km) and at 
a mountain (2770 m a.s.l). The observed results were compared with Monte Carlo calculations to 
find that an interaction model (Lund Fritiofl.6) used in an old neutrino flux calculation was not 
good enough for describing the observed values. In stead, we found that two other nuclear interac- 
tion models, Lund Fritiof7.02 and dpmjet3.03, gave much better agreement with the observations. 
Our data will serve for examining nuclear interaction models and for deriving a reliable absolute 
atmospheric neutrino flux in the GeV region. 



I. INTRODUCTION 

The discovery of evidence for neutrino oscillation by 
the Super Kamiokande group (l] is based on the compari- 
son of the observed atmospheric neutrino flux with calcu- 
lated values. Although the conclusion is so derived that 
it would not be upset by the uncertainty of the absolute 
flux value, it is desirable to obtain a reliable expected 
neutrino flux (under no oscillation assumption) for fur- 
ther detailed discussions. 

Two major sources of uncertainty in the atmospheric 
neutrino flux calculation are 1) the primary cosmic-ray 
spectrum and 2) the propagation of cosmic rays in the 
atmosphere, especially, modeling of the nuclear interac- 
tion. The absolute flux calculations so far made by vari- 
ous groups are expected to have uncertainty of ~ 30 %(2) . 

The primary proton and He spectra recently measured 
with magnet spectrometers by the BESS (|] and AMS|| 
groups agree very well and seem reliable. Therefore, we 
may take that the first problem mentioned above have 
now been almost settled at least up to 100 GeV/n. This 
means that if we have a reliable atmospheric cosmic-ray 
flux data, we may compare it with a calculation which 
uses such primaries and test the validity of nuclear inter- 
action models. 

For such an atmospheric cosmic-ray component, one 
may first raise the muon and actually some new observa- 
tions have been or being tried || q, [?]]. 

As a secondary cosmic-ray component, we focused on 
gamma-rays which are easy to measure with our dctcc- 
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tor. A good model should be able to explain muons 
and gamma-rays simultaneously. Muons arc important 
since they are directly coupled with neutrinos, but the 
flux is affected somehow by the structure of the atmo- 
sphere which is usually not well known. Compared to 
muons, the flux of gamma-rays is substantially lower but 
is almost insensitive to the atmospheric structure and 
depends only on the total thickness to the observation 
height. 

In 1998, we performed first gamma-ray observation 
with our detector at Mt. Norikura (2770m a.s.l) in Japan, 
and also made subsequent two successful observations at 
balloon altitudes (15 ~ 25 km) in 1999 and 2000. In the 
present paper, we report the final results of these obser- 
vations and consequences. 



II. THE DETECTOR 

For our observation, we upgraded the BETS (Balloon- 
born Electron Telescope with Scintillating fibers) detec- 
tor which had been developed for the observation of cos- 
mic primary electrons in the 10 GeV region. Its details 
before being upgraded for gamma-ray observation is in 
H and the electron observation result is in || . The basic 
performance was tested at CERN using electron, proton 
and pion beams of 10 to 200 GcV§, [UJ. Although this 
was undertaken before the upgrading, we can essentially 
use that calibration for the current observeions partly 
with a help of Monte Carlo simulations. 

Figure |l| shows a schematic structure of the main body 
of BETS. The calorimeter has 7.1 r.l lead thickness and 
the cross-section is 28 cm x 28 cm. The whole detector 
system is contained in a pressure vessel made of thin 
aluminum. 
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FIG. 1: Schematic illustration of the main body of the detec- 
tor. SI, S2 and S3 are 1 cm thick plastic scintillators used 
for trigger. Each fiber has 1mm diameter. Originally nu- 
clear emulsion plates were placed on the upper scifi's and also 
inserted between the upper thin lead plates for detailed inves- 
tigation of tracking capability of scifi. They are kept in the 
present system to have the same structure at the calibration 
time. The inlaid cascade shows charged particle tracks by a 
simulation for a 30 GeV incident proton. 



TABLE I: Basic characteristics of BETS 



(triple numbers in the table are for gamma-ray energy of 5, 
10, and 30 GeV, respectively) 



R.M.S energy resolution(%) 


21, 18, 15 (for ~ 15°) 


Sf2(cm 2 sr) 


243, 240,218 (at ~20 km) 


R.M.S angular resolution (deg) 


2.3, 1.3, 1.0 (for 6 ~ 15°) 


Total number of scifi's 


10080 


Weght including electronics (kg) 


230 


Cross-section of the main body 


28cm x 28cm 


Thickness (Pb radiation length) 


7.1 



R protons) at 10 GeV|2j| The basic characteristics of the 
detector are summarized in Table fl 
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FIG. 2: Image of cascade shower by a proton (120 GeV,left) 
and an electron(10 GeV, right) obtained at CERN. 

In Fig.||, we show examples of the CCD image of a cas- 
cade shower for a proton incident case and for an electron 
incident case. 

Figure || illustrates the yearly change of anti-counters. 
In 1998 (Mt.Norikura observation), the main change was 
limited to the upgrading of trigger logic. In 1999, we 
added 4 side anti-counters (each 15 cm x 36 cm x 1.5 cm 
plastic scintillator. Nine optical fibers containing wave 
length shifter are embedded in each scintillator and con- 
nected to a Hamamatu H6780 PMT. 




FIG. 3: Yearly change of the anti-counters. Left: 1998. No 
change from original BETS except for trigger logic. Middle: 
1999. 1.5 cm thick plastic scintillator side anti-counters were 
added. Right: 2000. The whole top view was covered by a 1 
cm thick plastic sintillator. 



The main feature of the BETS detector is that it is a 
tracking calorimeter; it contains a number of sheets con- 
sisting of 1 mm diameter scintillating fibers (scifi), many 
of which are sandwiched between lead plates. The to- 
tal number of scifi's arc 10080. The sheets are grouped 
into two types; one is to serve for x and the other for y 
position measurement. Each of them is fed to an im- 
age intensifier which in turn is connected to a CCD. 
Thus, the two CCD output gives us an x — y image of 
cascade shower development and enables us to discrimi- 
nate gamma-rays, electrons from other (mainly hadronic) 
background showers. The proton rejection power against 
electron is R ~ 2 x 10 3 (i.e, one misidentification among 



In 2000, we further added an anti-counter which covers 
the whole top view of the detector and also improved data 
acquisition speed. The top anti-counter is 38 cm x 38 
cm x 1 cm plastic scintillator. We also embedded optical 
fibers; 8 in the x and another 8 in the y direction, all of 
which were fed to an H6780. 

Although we could remove background showers with- 
out the anti-counters, inclined particles (mainly protons) 
entering from the gap between top scintillator (SI) and 
the main body degrades the desired gamma-ray event 
rate. The addition of the top anti-counter greatly helped 
improve this rate. 

We emphasize that detection of gamma-rays is easier 
for us than that of electrons, since, for gamma-rays, we 



3 



can utilize absence of incident charge. 



III. OBSERVATIONS 

Table |n] shows the summary of the observations. 

• Mt. Norikura observation. 

Our first gamma-ray observation was performed 
in 1998 at Mt. Norikura Observatory of Univ. of 
Tokyo, Japan (2770 m a.s.l, latitude 36.1°N, longi- 
tude 137.55°E, magnetic cutoff rigidity ~ 11.5 GV). 
The atmospheric pressure during the observation is 
shown in Fig.^. The average atmospheric depth is 
737 g/cm 2 . 
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FIG. 4: Pressure change during Mt. Norikura observation. 
The last pressure drop is due to a typhoon. The average 
pressure is 723 hP (737 g/cm 2 ). 



• Balloon flight 

We had two similar balloon Alights in 1999 and 
2000. Since the main outcome of the data is from 
the latter, we briefly describe it. A balloon of 
43 x 10 3 m 3 was launched at 6:30 am, 5th June, 2000 
from the Sanriku balloon center of the Institute of 
Space and Astronautical Science, Japan (latitude 
39.2°N, longitude 141. 8°E, magnetic cutoff rigidity 
~ 8.9 GV) and recovered with the help of the heli- 
copter, at 17:59 on the sea not far from the center. 
The flight curve shown in Fig.|5| confirms that we 
have good level flights at 4 different heights. 

As compared to the 1999 flight, this flight real- 
ized a smaller dead time and higher ratio of desired 
gamma-ray events. 
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FIG. 5: Flight curve of the 2000 observation. Pressure (up- 
per) and altitude (lower) as a function of time. Each arrow 
shows the level flight region. The pressure change at around 
15.3 km is rather rapid but the gamma-ray intensity is almost 
constant there and the change can be neglected. 



A. Event trigger 

The basic event trigger condition is created by signals 
from the three plastic scintillators (SI, S2 and S3). We 
show the discrimination level in terms of the minimum 
ionizing particle number which is defined by the peak of 
the energy loss distribution of cosmic-ray muons passing 
both SI and S3 with inclination less than 30 degrees. 

We prepare a multi-trigger system by which event trig- 
ger with different conditions is possible at the same time. 
The major two trigger modes are the g-low and g-high. 
The g-low is responsible for low energy gamma-rays and 
all anti-counters, when available, are used as veto coun- 
ters. Its condition is listed in Table O. High energy 
gamma-rays normally produce a lot of back splash par- 
ticles which hit SI and/or anti-counters, and thus the 
g-low trigger is suppressed. In such a case, i.e, if we have 
a large S3 signal, anti-counter veto is invalidated and the 
SI threshold is relaxed (The g-high condition is Sl< 3.0, 
S2> 5.0 and S3> 8.1). 

The branch even point of the g-low and g-high mode 
efficiency is at ^30 GcV. Since we deal with gamma-rays 
mostly below 30 GeV, and also to avoid complexity, we 
present results only by the g-low mode. 
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TABLE II: Summary of three observations 



Observation 


Mt.Norikura(1998) 


Balloon (Sanriku, 1999) 


Balloon (Sanriku, 2000) 


Period 


Aug.31~Sep.18 




Sep.2, 6:55- 


-17:17 






Jun.5, ( 


>:30~17:59 


Altitude(km) 


2.77 


15.3 


18.5 


21.2 


24.7 


32.3 


15.3 


18.3 


21.4 


25.1 


Depth(g/cm 2 ) 


737 


126 


74.8 


48.9 


28.0 


9.5 


128 


73 


45.7 


25.3 


Obs. hour (s) 


1.33 x 10 6 


1260 


1560 


2100 


4878 


3120 


1560 


2160 


4320 


2320 


Live time (s) 


9.8 x 10 5 


504 


450 


414 


852 


498 


752 


928 


1805 


789 


Live time (%) 


74.0 


40.0 


28.8 


19.7 


17.5 


16.0 


48.2 


43.0 


42.6 


44.2 


Triggered events 


1.8 x 10 


9513 


11288 


13361 


30439 


16741 


18808 


25795 


46675 


17436 


7 events 


4.7 x 10 4 


700 


650 


611 


848 


345 


1300 


1485 


2299 


740 


(%) 


2.5 


7.3 


5.7 


4.6 


2.8 


2.0 


6.9 


5.8 


4.9 


4.2 


g-low trigger 


SI < 0.5 






Sl< 0.5 






Sl< 0.47 




condition (in mip). 


S2 > 2.3 






S2> 1.5 






S2> 1.59 






S3 > 1.7 






S3> 3.0 






S3> 3.18 





IV. ANALYSIS 



Event selection 



Among the triggered events, we selected gamma-ray 
candidates by imposing the following conditions: 



analysed events together with the result of CERN data. 
Hadrons make a distribution with a peak at around 0.5. 
We see that the contribution of hadrons in our observa- 
tion is negligible. 



B. Energy Determination 
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FIG. 6: (left)Energy concentration distribution at 21.4 km. 
(right)the same by electrons at CERN 



1. The estimated shower axis passes SI and S3. The 
axis position in S3 must be at least 2 cm apart from 
the edge of S3. 

2. The estimated shower axis has a zenith angle less 
than 30 degrees. 

3. The energy concentration (see below) must be 
greater than 0.7. 

According to a simulation, only neutrons could be a 
background against gamma-rays and the 3rd conditions 
above reduces the neutron contribution to a negligible 
level (< 1%). 

The energy concentration is defined as the fraction 
of scintillating fiber light intensity within 5 mm from 
the shower axis. Figure ^| shows the concentration of 



The energy calibration was performed in 1996 at 
CERN using electrons with energy 10 - 200 GcV[|, |o|. 
There is no direct calibration for gamma-rays, but, for 
the present detector thickness and energy range, a M.C 
simulation tells us that the calibration in 1996 can be 
used for gamma-rays, too[^4). Therefore, for the 1998 
and 1999 observations, energy is obtained as a function 
of the S3 output and zenith angle using the CERN cali- 
bration. 

In 2000, we made some change in the electronics so 
the CERN calibration could not be used directly. The 
effect by the change was absorbed by a M.C simulation 
of which the validity was verified by examining the 1998 
and 1999 data. We used the sum of S2 and S3 outputs 
below 20 GeV since the energy resolution was found to be 
better than using S3 only. Figure |?] shows r.m.s energy 
resolution. 



C. Correction of the gamma-ray intensity 

The gamma-ray vertical flux is obtained from the raw 
dN/dE by dividing it by the live time of the detector 
and the effective 5f2 (area x solid angle). The latter 
is obtained by a simulation ]ll[. It is dependent on the 
observation hight and energy. A typical value at 10 GeV 
is 240 cm 2 sr (see Table|). The energy spectrum is further 
corrected by the following factors which are not taken 
into account in the SQ calculation. 
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FIG. 7: R.m.s energy resolution. The resolution by S2+S3 or 
S3 only is shown. Different symbols indicate different incident 
angles. We used S2+S3 below 20 GeV for the year 2000 data. 



Systematic bias in our estimation of the shower 
axis. We underestimate the zenith angle systemat- 
ically and it leads to overestimation of the intensity 
about 4% for the balloon and 1.8 % for Mt.Norikura 
observations. 

Multiple incidence of particles. A gamma-ray is 
sometimes accompanied by other charged particles 
and they enter the detector simultaneously (within 
1 ns time difference in 99.9 % cases). They are a 
family of particles generated by one and the same 
primary particle ||f|. The charged particles fire the 
anti-counter and the g-low trigger is inhibited. 

In some case, multiple gamma-rays enter the de- 
tector simultaneously. The rate is smaller than the 
charged particle case. However, this is judged as 
a hadronic shower in most of cases. The multiple 
incidence leads to the underestimation of gamma- 
ray intensity. The portion of multiple incidence is 
shown in Fig.^J (upper). 

Finite energy resolution. The rapidly falling en- 
ergy spectrum leads to the spillover effect. This 
normally leads to the overestimation of flux (FigJ|, 
lower) . 



RESULTS AND COMPARISON WITH 
CALCULATIONS 



The flux values are summarized in Table II]. We put 
only the statistical errors in the flux values, since system- 
atic errors coming from the uncertainty of the Sf2 calcula- 
tion, various cuts and flux corrections are expected to be 
order of a few percent and much smaller than the present 
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FIG. 8: (upper) Multiple incidence rate, (lower) Correction 
factor for year 2000 due to spillover. The flux must be low- 
ered. For Norikura, the factor below 20 GeV is larger by 1~ 3 
%. 



statistical errors. 

The gamma-ray energy spectra thus obtained at bal- 
loon altitudes are shown in Fig.|| together with the 
expected ones calculated by the Cosmos simulation 
Except for 32.3 km altitude, we can disregard 
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the small difference of the observation depths and we 
combine two flight data with statistical weight, although 
the main contribution is from the flight in 2000. 

In the simulation calculation, we employed 3 different 
nuclear interaction models: 1) fritiof 1 . 6 [p6| used in 
the HKKM calculation jl^], which was widely used for 
comparison with the Kamioka data, 2)fritiof7.02 [|l5|] |27[| 
and 3) dpmjet3.03jl6|. As the primary cosmic ray, we 
used the BESS result on protons and He. The CNO com- 
ponent is also considered Q. Besides these we included 
electron and positron data by AMS]lq|. Their data in 
the 10 GeV region is consistent with the HEAT[jl9| and 
BETS H data. Bremstrahlung gamma-rays from the pri- 
mary electrons could contribute order of ~ 10 % at very 
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TABLE III: Summary of flux values 



height (km) 


15.3 


18.3 


21.4 


25.1 


32.3 








Energy (GeV) & flux(No./m 2 - 


s-sr-GeV) 






5.48 


2.42 ± 0.37 


5.48 


2.11 ± 0.39 


5.47 


2.11 ± 0.24 


5.47 


1.58 ± 0.25 


5.47 


0.49 ± 0.14 


6.47 


1.18 ± 0.27 


6.47 


1.10 ± 0.24 


6.47 


1.35 ± 0.21 


6.47 


0.82 ± 0.18 


6.57 


0.19 ± 0.09 


7.47 


0.89 ± 0.24 


7.47 


0.79 ± 0.21 


7.47 


0.82 ± 0.16 


7.47 


0.66 ± 0.16 


7.47 


0.24 ± 0.10 


8.48 


0.37 ± 0.15 


8.48 


0.92 ± 0.20 


8.48 


0.51 ± 0.13 


8.48 


0.49 ± 0.14 


8.48 


0.16 ± 0.08 


9.48 


0.54 ± 0.17 


9.85 


0.46 ± 0.11 


9.48 


0.50 ± 0.12 


9.48 


0.36 ± 0.12 


9.48 


0.16 ± 0.08 


10.5 


0.17 ± 0.10 


11.5 


0.35 ± 0.12 


10.5 


0.41 ± 0.09 


10.5 


0.34 ± 0.12 


12.3 


0.13 ± 0.037 


12.1 


0.28 ± 0.09 


14.0 


0.24 ± 0.06 


11.8 


0.23 ± 0.069 


12.2 


0.21 ± 0.054 


17.0 


0.032 ± 0.018 


14.0 


0.17 ± 0.05 


18.3 


0.072 ± 0.030 


14.0 


0.16 ± 0.030 


14.0 


0.076 ± 0.03 


21.7 


0.022± 0.015 


18.5 


0.12 ± 0.04 


26.8 


0.040 ± 0.017 


18.4 


0.086 ± 0.023 


17.8 


0.078 ± 0.029 






25.5 


0.06 ± 0.02 






27.1 


0.026 ± 0.009 


21.7 
26.8 
36.0 


0.064 ± 0.026 
0.024 ± 0.012 
0.012 ± 0.008 







TABLE IV: Flux values at Mt. Norikura 



E(GeV) 


Flux (10" 4 /m 2 -s-sr-GeV) 


5.48 


274 ± 13 


6.47 


183 ± 11 


7.47 


133 ± 9 


8.47 


87.8 ± 7.5 


9.47 


86.5 ± 7.5 


10.5 


54.1 ± 5.9 


11.5 


46.6 ± 5.5 


12.5 


38.3 ± 5.0 


13.5 


32.6 ± 4.6 


14.5 


24.2 ± 4.0 


15.5 


25.7 ± 4.1 


17.0 


11.9 ± 2.0 


19.0 


15.3 ± 2.3 


21.0 


13.1 ± 2.1 


23.0 


5.80 ± 1.4 


26.0 


5.31 ± 0.95 


30.0 


3.00 ± 0.72 


34.0 


2.30 ± 0.64 


38.0 


1.07 ± 0.44 


45.0 


1.45 ± 0.32 


55.0 


0.52 ± 0.20 


65.0 


0.22 ± 0.13 


75.0 


0.30 ± 0.15 


85.0 


0.15 ± 0.10 



high altitudes. 

At balloon altitudes, the two models, fritiof7.02 and 
dpmjct3.03, give almost the same results which are close 
to the observed data, while fritofl.6 gives clearly smaller 



fluxes than the observation. 

Figure [l(] shows the result from the observation at 
Mt. Norikura. It should be noted that the flux by 
fritiofl.6 becomes higher than the ones by the other mod- 
els at this altitude. 

From these figures, we see fritiof7.02 and dpmjet3.03 
give rapider increase and faster attenuation of intensity 
than fritiofl.6; the tendency is very consistent with the 
observed data. The transition curve of the flux integrated 
over 6 GeV shown in Fig.jll] clearly demonstrates this 
feature. 
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FIG. 9: Gamma-ray spectra at 5 balloon heights are compared 
with 3 different models. The vertical axis is FluxxS 2 . Ex- 
cept for 1999 data at 32.3 km, 1999 and 2000 flights data are 
combined. From top to bottom, at 25.1, 21.4, 18.3, 15.3 and 
32.3 km. The spectra expected from three interaction models 
are drawn by solid (dpmjet3.03), dash (fritiof7.02) and dotted 
(fritiofl.6) lines. 
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FIG. 10: Gamma-ray spectrum at Mt. Norikura (2.77 km 
a.s.l). The vertical axis is Fluxxi? 2 . Our data is at < 100 
GeV. Data above 300 GeV is from emulsion chamber experi- 
ments. For the latter, see Sec. VI 
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FIG. 11: The altitude variation of the flux integrated over 
6 GeV. The dpmjet3.03 and fritiof7.02 give almost the same 
feature consistent with the observation while the deviation of 
fritiofl.6 from the data is obvious. 



VI. DISCUSSIONS 

A. Comparison with other data 

We found Fritiof7.02 and dpmjet3.03 give good agree- 
ment with the observed gamma-ray data at around 10 
GeV. We briefly see whether these models can interpret 
other observations. More detailed inspection will be done 
elsewhere. 

• Muon data by the BESS group at Mt.NorikuraQ. 

Recently, the BESS group reported detailed muon 
spectrum over several hundred McV/c. In their 



paper, calculations by dpmjet3.03 and fritiofl.6 are 
compared with the data; agreement by dpmjet3.03 
is quit good at least above GeV where Fritiof7.02 
also gives more or less the same flux. On the other 
hand, fritiofl.6 shows too high flux. These features 
are consisten with our present analysis. 

• Higher energy gamma-ray data by emulsion cham- 
ber. 

In Fig. nfl , we inlaid an emulsion chamber 
datafl20|| p8| at Mt. Norikura. Our data seems to 
be smoothly connected to their data as the two in- 
teraction models (Fritiof7.02 and dpmjet3.03) pre- 
dict. Since the emulsion chamber data extends to 
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the TeV region and the primary particle energy re- 
sponsible for such high energy gamma-rays is much 
higher than 100 GeV where we have no accurate in- 
formation comparable to the AMS and BESS data, 
it would be premature to draw a definite conclusion 
on the primary and interaction model separately. 
However, the fact that smooth extrapolation of the 
primary spectra as shown in Table |V] and the in- 
teraction model, dpmjet3.03 or fritiof7.02, give a 
consistent result with the data, seems to indicate 
that such combination would provide a good esti- 
mate on other components at 3> 10 GeV. 



TABLE V: Primary flux assumed in the simulation above 100 
GeV/n 

(E in kinetic energy per nucleon (GeV), flux in /m 2 -s-sr-GeV) 



Proton 


Helium 


CNO 


E 


flux 


E flux 


E 


flux 


92.6 


0.593E-01 


79.4 0.549E-02 


100. 


9.0E-5 


108 


0.388E-01 


100. 3.0E-3 


400. 


1.8E-6 


126 


0.276E-01 


200. 5.0E-4 


2.0E3 


3.5E-8 


147 


0.179E-01 


400. 7.0E-5 


2.0E4 


9.3E-11 


171 


0.124E-01 


2.0E3 9.98E-7 


2.0E5 


2.3E-13 


200 


0.836E-02 


2.0E4 2.5E-9 


14.0E5 


1.3E-15 


1100 


8.29E-5 


2.0E5 3.97E-12 


3.0E6 


1.7E-16 


1.1E4 


1.47E-7 


4.0E5 6.1E-13 


3.0E7 


2.0E-19 


1.1E5 


2.8E-10 


8.0E5 7.0E-14 


3.0E8 


2.2E-22 


2.2E5 3.7E-11 


8.0E6 8.7E-17 






4.4E5 


5.0E-12 


8.0E8 5.3E-23 






4.4E8 


2.8E-21 










B. 


The ^-distributions 





The two models, fritiof7.02 and dpmjet3.03, give al- 
most the same results in the present comparison. How- 
ever, if we look into the x-distribution of the particle pro- 
duction, we note some difference, especially in the proton 
x-distribution. We define the x as the kinetic energy ra- 
tio of the incoming proton and a secondary particle in 
the laboratory frame. The x distribution for pAir col- 
lisions at incident proton energy of 40 GeV is presented 
for photons (from n° plus r] decay) and protons in Figjl^. 
Difference of the three models seen in the photon distri- 
bution is quite similar to the one for charged pions. The 
x region most effective to atmospheric gamma-ray flux is 
around 0.2~0.3 where the difference is not so large but 
fritiof7.02 and dpmjct3.03 have higher gamma-ray yield 
than fritiofl.6. 

On the other hand, the proton x distribution has larger 
difference among the three models (we note, however, 
the difference may be exaggerated than the photon case 
due to the scale difference). It is interesting to sec 
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FIG. 12: The x-distribution of photons from 7r° plus r\ decay 
(upper) and protons (lower) for pAir collisions at 40 GeV. 
The three model results are shown. 



that, in spite of these large differences, the final flux is 
not so much different each other. Our gamma-ray data 
prefers to rather more inelastic feature of collisions than 
fritiofl.6, i.e rapider increase and faster attenuation of 
the flux. 

We should compare the distribution with accelerator 
data; however, there is meager stuff appropriate for our 
purpose. One such comparison has been done in a recent 
review paper Q for pAir collisions at 24 GeV/c incident 
momentum. The charged pion distribution by fritiofl.6 
and dpmjet3.03 well fit to some scattered data which 
prevents to tell the superiority of the two. As to the 
proton distribution, among the three models, fritiofl.6 is 
rather close to the data but deviation from the data is 
much larger than the pion case. 

The proton ^-distribution would strongly affect the at- 
mospheric proton spectrum. We calculated proton flux at 
Mt.Norikura to find a flux relation such that fritiofl.6 > 
fritiof7.02 > dpmjet3.03 as expected naturally from the 
^-distributions. The maximum difference is factor ~ 2.5 
in the energy region of 0.3 to 3 GeV. The BESS group 
has measured the proton spectrum at Mt. Norikura in 
the same energy region. Their result expected to come 
soon[|l) will help select a better model for the proton x 
distribution. 
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VII. SUMMARY 

• We have made successful observation of at- 
mospheric gamma-rays at around 10 GcV at 
Mt.Norikura (2.77 km a.s.l) and at balloon alti- 
tudes (15 - 25 km). 

• The observed gamma-ray fluxes are compared with 
calculations by three interaction models; it is found 
that fritiofl.6 employed by the HKKM calcula- 
tion fll4[ , which was used in comparison with the 
Kamioka data, is not a very good model. 

• Other two models (fritiof7.02 and dpmjet3.03) give 
better results consistent with the data, which shows 
rapider increase and faster attenuation of the flux 
than fritiofl.6 predicts. 

• Our data has complementary feature to muon data 
and will serve for checking nuclear interaction mod- 
els used in atmospheric neutrino calculations. 
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